System and method for assessing wound

ABSTRACT

The wound assessing method and system of the present teachings provide a convenient, quantitative mechanism for diabetic foot ulcer assessment.

CROSS REFERENCE TO RELATED APPLICATIONS

This application claims priority to and benefit of U.S. ProvisionalApplication No. 61/897,559, entitled SYSTEM AND METHOD FOR ASSESSINGWOUND, filed on Oct. 30, 2013, and of U.S. Provisional Application No.61/898,907, entitled SYSTEM AND METHOD FOR ASSESSING WOUND, filed onNov. 1, 2013, both of which are incorporated by reference herein intheir entirety and for all purposes.

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT

This invention was made with government support under Grant No.IIS-1065298, awarded by the National Science Foundation (NSF). Thefederal government may have certain rights in the invention.

BACKGROUND

The present teachings relate to a system and a method for assessingchronic wound. More particularly, the present teachings relate to asystem and a method for assessing wound for patients with, for example,type 2 diabetes and diabetic foot ulcers. One way to assess wound is touse a specialized camera to capture the wound image, then calculates thewound area and organizes wound images from different patients and storesimages in a central location. Another way to assess wound is to use amobile wound analyzer (MOWA), which is an Android-based software,intended for smart phones and tablets, for analysis of wound images. Thewound boundary needs to be traced manually after which the softwarecalculates the wound area and performs color analysis within the woundboundaries.

The conventional art does not address the problem of capturing footimages when the patients with diabetes have limited mobility. Inaddition, the prior art device is very costly and not affordable forindividual patients to own, apart from MOWA, which, however, is designedfor clinicians. Further, the prior art is not designed for joint use byboth the patient and his/her doctor, through automatic upload of raw andanalyzed wound images to cloud storage for easy access by the physician.Accordingly, there is a need to develop new system and method forassessing wound that overcome the above drawbacks in the prior art.

BRIEF SUMMARY

The present teachings provide patients with diabetes and chronic footulcers an easy-to-use and affordable tool to monitor the healing oftheir foot ulcers via a healing score; at the same time, the patient'sphysician can review the wound image data to determine whetherintervention is warranted. The system is also applicable for patientswith venous leg ulcers.

In accordance with one aspect, the present teachings provide a methodfor assessing wound. In one or more embodiments, the method of theseteachings includes capturing an image of a body part including the woundarea, analyzing the image to extract a boundary of the wound area,performing color segmentation within the boundary, wherein the woundarea is divided into a plurality of segments, each segment beingassociated with a color indicating a healing condition of the segmentand evaluating the wound area.

In accordance with another aspect, the present teachings provide asystem for assessing wound. In one or more embodiments, the system ofthese teachings includes an image acquisition component configured tocapture an image of a body part including a wound area, an imageanalysis module configured to extract a boundary of the wound area; animage segmentation module configured to perform color segmentationwithin the boundary of the wound area, wherein the wound area is dividedinto a plurality of segments, each segment being associated with a colorindicating a healing condition of the segment and a wound evaluationmodule configured to evaluate the wound area.

A number of other embodiments of the method and system of theseteachings are presented herein below.

For a better understanding of the present teachings, together with otherand further needs thereof, reference is made to the accompanyingdrawings and detailed description and its scope will be pointed out inthe appended claims.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 represents a schematic block diagram representation of oneembodiment of the system of these teachings;

FIG. 2a is a flowchart representation of the mean shift basedsegmentation algorithm as used in these teachings;

FIG. 2b is a flowchart representation of the mean shift based method forwound boundary determination as used in these teachings;

FIGS. 3a-3d are results of different operations in the wound boundarydetermination of these teachings;

FIGS. 4a-4c are results in the calculation of wound location of theseteachings;

FIG. 5 is a flowchart representation of the machine learning basedmethod for wound recognition of these teachings;

FIG. 6 shows exemplary results of a wound image obtained by a majorityvote scheme as used in these teachings;

FIG. 7 shows a comparison of original images, images obtained by themajority vote scheme and wound recognition using the machine learningmethods of these teachings;

FIG. 8 shows the conventional confusion matrix that is used in theseteachings;

FIG. 9a is a flowchart representation of the K-means algorithm used inthese teachings;

FIG. 9b is a flowchart representation of Color segmentation method ofthese teachings;

FIG. 9c is a flowchart representation of algorithm flow of the regionfusion step of these teachings;

FIGS. 10a-10i show images of wound areas and all results of two colorsegmentation methods used in these teachings;

FIG. 11 is a schematic representation of a component of one embodimentof the system of these teachings;

FIGS. 12a-12c are graphical representations of a component of oneembodiment of the system of these teachings;

FIG. 13 is a schematic representation of front surface mirrors as usedin one component of the system of these teachings;

FIG. 14 is a schematic block diagram representation of anotherembodiment of the system of these teachings;

FIG. 15 is a schematic block diagram presentation of another componentof a further embodiment of the system of these teachings; and

FIG. 16 is a flowchart representation of the use of a further embodimentof the system of these teachings.

DETAILED DESCRIPTION

The following detailed description is not to be taken in a limitingsense, but is made merely for the purpose of illustrating the generalprinciples of these teachings, since the scope of these teachings isbest defined by the appended claims. Although the teachings have beendescribed with respect to various embodiments, it should be realizedthese teachings are also capable of a wide variety of further and otherembodiments within the spirit and scope of the appended claims.

As used herein, the singular forms “a,” “an,” and “the” include theplural reference unless the context clearly dictates otherwise.

Except where otherwise indicated, all numbers expressing quantities ofingredients, reaction conditions, and so forth used in the specificationand claims are to be understood as being modified in all instances bythe term “about.”

In the following, the term “handheld mobile communication device,” asused herein, refers to a device capable of being handheld and ofexecuting applications, and which is portable. In one instance, themobile communication device has one or more processors and memorycapability. Examples of mobile communication devices, these teachingsnot being limited to only these examples, include smart mobile phones,digital personal assistants, etc.

The present teachings relate to a wound image analysis system, which maybe implemented as hardware and/or software. In one embodiment, the woundimage analysis system of the present teachings is designed to operate ona handheld mobile communication device, such as a smart phone. The woundimage analysis system may be used in private homes or elder carefacilities by the patient him/herself, or in collaboration with acaregiver, with the relevant image data automatically uploaded to securecloud storage, to be accessible for perusal by the patient's doctorand/or clinicians in the patient's wound clinic. An alert system cannotify the patient's doctor if wound data exceeds some preset bounds. Inanother embodiment, the wound image analysis system of the presentteachings may operate in a wound clinic and cover several patients. Inthis embodiment, a smart phone is use collaboratively with a laptop(i.e., a smart phone-laptop collaborative system).

In one embodiment, the wound image analysis method of the presentteachings includes the following main parts: (i) image preprocessing,(ii) method for determining the wound boundary, (iii) method for colorimage segmentation, (iv) method for computing the healing score. Inother embodiments, the system of these teachings component configured todetermine the wound boundary, configured to perform color imagesegmentation and component configured to assess the wound area. Otherembodiments of the system of these teachings also include an imagecapture box to aid the patient and/or his/her caregiver in capturingimages of the foot ulcer under controlled distance and light conditions,and cloud storage and clinical access solution. Each of these componentswill be described briefly below, with additional details given in theattached documents. While each system component is essential for thefunctionality of the system, not all components are necessary to operatethe wound image analysis system.

(i) Image pre-processing. A JPEG image captured by a smart phone isconverted into an RUB bitmap image. An image noise reduction filter isapplied to down-sample the image for faster processing.

(ii) Component configured to determine the wound boundary. The woundboundary detection method is based on the mean shift segmentation of thewound image. The method first detects the outline of the foot and thenwithin the boundary of the foot locates the outline of the wound. A moreaccurate method may be used for wound boundary detection based on skillsand insight by experienced wound clinicians. For this purpose, machinelearning methods, such as the Support Vector Machine, may be used totrain the wound analysis system to learn about the essential featuresabout the wound.

(iii) Component configured for color image segmentation. The colorsegmentation method is instrumental in determining the healing state ofthe wound where red indicates healing, yellow indicates inflamed, andblack indicates necrotic.

(iv) Component configured to compute a healing score. The Healing Scoreis an important element of communicating in a simple fashion the healingstatus of the patient's wound. The Healing Score is a weighted sum offactors, such as: wound area; weekly change in wound area; woundtexture; relative size and shapes of the healing, inflamed and necroticregions within the wound boundary, and possibly the skin color aroundthe wound. The weighing factors are determined from expert clinicalinput.

(v) Image capture box. The image capture box is a device that allows apatient, possibly with the aid of his/her caregiver, to both visuallyobserve the appearance of a wound on the sole of the foot as well ascapture an image of the wound for storage and analysis. It is a compactbox, where the patient's foot can rest comfortably on a 45′ angledsurface next to the smart phone holder. The angled surface can readilybe flipped to accommodate right foot as well as left foot. The boxcontains two front surface mirrors and warm white LED lighting.

(vi) Cloud storage and clinical access solution. The cloud storage andclinical access solution automatically uploads relevant wound data tothe cloud (e.g., network accessible storage) from the smart phone,either utilizing Wi-Fi (802.11), 3G, or other wireless network. Relevantdata comprises wound image data, which is automatically uploaded inencrypted form to secure cloud storage, to be accessible for perusal bythe patient's doctor. An alert system can alert the doctor if wound dataexceeds some preset bounds.

In another embodiment, the wound image analysis system operates in awound clinic and covers several patients. In this embodiment, a handheldmobile communication device-Computing Component collaborative system isused, in which a captured image is automatically transferred to acomputing component, in one instance, the transfer occurs through apeer-to-peer based Wi-Fi system or Local Area Network, using a wired orwireless router.

Moreover, in a further embodiment, instead of using the smart phonecamera as the image acquisition device, the wound image analysis systemcan use a compact hyperspectral camera integrated into the image capturebox. In one instance, three types of LED illumination are integratedinto the image capture box: infrared (IR) LED illumination; visiblelight illumination, using the already built-in warm white LEDillumination; and ultraviolet (UV) LED illumination. This allows thewound to be imaged by three distinct wavelength bands, with theexpectation of revealing much better diagnostic information about thewound. In one instance, The hyperspectral camera includes directcommunication capability, such as but not limited to by which thecaptured images are transmitted to a device, such as a handheld mobilecommunication device or a computing device, for processing and cloudupload.

In accordance with one aspect, the present teachings provide a methodfor assessing wound. In one or more embodiments, the method of theseteachings includes capturing an image of a body part including a woundarea, analyzing the image to extract a boundary of the wound area,performing color segmentation within the boundary, wherein the woundarea is divided into a plurality of segments, each segment beingassociated with a color indicating a healing condition of the segmentand evaluating the wound area.

In one or more embodiments, the system of these teachings includes animage acquisition component configured to capture an image of a bodypart including a wound area, an image analysis module configured toextract a boundary of the wound area; an image segmentation moduleconfigured to perform color segmentation within the boundary of thewound area, wherein the wound area is divided into a plurality ofsegments, each segment being associated with a color indicating ahealing condition of the segment and a wound evaluation moduleconfigured to evaluate the wound area.

One embodiment of the system of these teachings is shown in FIG. 1.Referring to FIG. 1, in the embodiment shown there in, an image capturecomponent 15 captures the image, the image is preprocessed and providedto an image analysis component 25 that is configured to extract aboundary of the area of the wound. An image segmentation component 35 isconfigured to perform color segmentation on the image within the woundboundary. In the color segmentation the area of the image within the oneboundary is divided into a number of segments, each segment beingassociated with the color that indicates a healing condition of thesegment. A wound evaluation component 45 receives the information fromthe image segmentation component 35 and provides an analysis of thewound healing trend.

In one instance, after the wound image is captured, the JPEG file pathof this image is added into a wound image database. This compressedimage file, which cannot be processed directly with the main imageprocessing algorithms, therefore needs to be decompressed into a 24-bitbitmap file based on the standard RGB color model. In one instance, thebuilt-in APIs of the smartphone platform to accomplish the JPEGcompression and decompression task. The “image quality” parameter isused to control the JPEG compression rate. In one embodiment, setting“image quality” to 80 was shown empirically to provide the desirablebalance between quality and storage space. In that embodiment, for anefficient implementation on the smartphone alone, no method was used tofurther remove the artifacts introduced by JPEG lossy compression.

In one instance, in the Image preprocessing step, the high resolutionbitmap image is first down-sampled to speed up the subsequent imageanalysis and to eliminate excessive details that may complicate thewound image segmentation. In one instance, the original image (pixeldimensions 3264×2448) is down-sampled by a factor 4 in both thehorizontal and vertical directions to pixel dimensions of 816×612, whichhas proven to provide a good balance between the wound resolution andthe processing efficiency. Afterwards, the images is smoothed to removenoise (assumed mainly to be Gaussian noise produced by image acquisitionprocess) by using the Gaussian blur method whose standard deviationσ=0.5 was empirically judged to be substantially optimal based onmultiple experiments.

In one or more instances, in the method of these teachings, analyzingthe image includes performing mean shift segmentation and objectrecognition and, in the system of these teachings, the image analysiscomponent is configured to perform mean shift segmentation and objectrecognition. In the mean shift based image segmentation and region mergeoperations, the wound boundary determination task doesn't rely on anyclinical inputs. The Foot outline detection is accomplished by findingthe largest connected component in the segmented image. Afterwards, awound boundary determination was carried out based on the smart analysisof previous foot outline detection result. This solution is veryefficient and easy to be implemented on the handheld mobilecommunication device platform. However, the toe-amputation status has tobe recorded as part of patients' data.

Basic Mean Shift Method

Many non-parametric clustering methods can be separated into two parts:hierarchical clustering and density estimation. Hierarchical clusteringis a method of cluster analysis, which seeks to build a hierarchy ofclusters. Strategies for hierarchical clustering generally fall into twotypes including 1) agglomerative: this a “bottom up” approach in whicheach observation starts in its own cluster and pairs of clusters aremerged as one moves up the hierarchy; 2) divisive: this is a “top down”approach in which all observations start in one cluster and splits areperformed recursively as one moves down the hierarchy. On the otherhand, the concept of the density estimation-based non-parametricclustering method is that the feature space can be considered as theexperiential probability density function of the represented parameter.The mean shift algorithm can be classified as density estimation. Itadequately analyzes feature space to cluster them and can providereliable solutions for many vision tasks.

In general, the mean shift algorithm models the feature vectorsassociated with each pixel (e.g., color and position in the image grid)as samples from an unknown probability density function ƒ(x) and thentry to find clusters (dense areas) in this distribution. The key to meanshift is a technique for efficiently finding peaks in thishigh-dimensional data distribution (In these teachings, there will be 5dimension including 3 color range dimension and 2 spatial dimension)without ever computing the complete function explicitly. The problem issimplified to how to find the local maxima (peaks or modes) in anunknown density distribution. Let us take a look at the kernel densityestimation definition at first. Given n data points x_(i), i=1, . . . ,n in the d-dimensional space R^(d), the multivariate kernel densityestimator with kernel K(x) is shown as below (see D. Comdniciu, P. Meer,Mean Shift: A Robust Approach Toward Feature Space Analysis, IEEE Tran.on Pattern Analysis and Machine Intelligence, vol 24 (5), May 2002, pp.603-619, which is incorporated by reference herein is entirety and forall purposes).

$\begin{matrix}{{f(x)} = {\frac{1}{{nh}^{d}}{\sum\limits_{i = 1}^{n}\;{K\left( \frac{x - x_{i}}{h} \right)}}}} & {{Eq}.\mspace{14mu} 3.1}\end{matrix}$where h is one bandwidth parameter satisfying h>0 and K is the radiallysymmetric kernels satisfyingK(x)=c _(k,d) k(∥x∥ ²)  Eq. 3.2where c_(k,d) is a normalization constant which makes K(x) integrate toone. The function k(x) is the profile of the kernel, defined only forx≥0. After applying the profile notation, the density estimator in Eq.3.1 can be written as below [32].

$\begin{matrix}{{f_{h,K}(x)} = {\frac{c_{k,d}}{{nh}^{d}}{\sum\limits_{i = 1}^{n}\;{k\left( {{\frac{x - x_{i}}{h}}}^{2} \right)}}}} & {{Eq}.\mspace{14mu} 3.3}\end{matrix}$

In mean shift algorithm, a variant of what is known in the optimizationliterature is used as multiple restart gradient descent. Starting atsome guess for a local maxima y_(k), which can be a random input datapoint x_(i), mean shift computes the density estimate ƒ(x) at y_(k) andtake a uphill step in the gradient descent direction. The gradient ofƒ(x) is given by

$\begin{matrix}{{\nabla{f(x)}} = {\frac{2c_{k,d}}{{nh}^{d + 2}}{\sum\limits_{i = 1}^{n}\;{\left( {x_{i} - x} \right){g\left( \frac{{{{x_{i} - x}}}^{2}}{h^{2}} \right)}}}}} & {{Eq}.\mspace{14mu} 3.4}\end{matrix}$where g(r)=−k′(r) and n is the number of neighbors taken into account inthe 5 dimension sample domain. In one instance, the Epanechinikov kernelshown as Eq. 3.2 is used, which makes the derivative of this kernel is aunit sphere.If the Eq. 3.4 is rewritten as the following

$\begin{matrix}{{\nabla{f(x)}} = {{\frac{2c_{k,d}}{{nh}^{d + 2}}\left\lbrack {\sum\limits_{i = 1}^{n}\;{g\left( {{\frac{x - x_{i}}{h}}}^{2} \right)}} \right\rbrack} \times {m(x)}}} & {{Eq}.\mspace{14mu} 3.5} \\{{m(x)} = \frac{\sum\limits_{i = 1}^{n}\;{x_{i}{g\left( {{\frac{x - x_{i}}{h}}}^{2} \right)}}}{\sum\limits_{i = 1}^{n}{g\left( {{\frac{x - x_{i}}{h}}}^{2} \right)}}} & {{Eq}.\mspace{14mu} 3.6}\end{matrix}$

The vector in Eq. 3.6 is called the mean shift vector, since it is thedifference between the weighted mean of the neighbors x_(i) around x andthe current value x. In the mean-shift procedure, the current estimateof the mode y_(k) at iteration k is replaced by its locally weightedmean as shown below:y _(k) =y _(k) +m(y _(k))  Eq. 3.7

This iterative update of the local maxima estimation will be continueduntil the convergence condition is met. In one instance, the convergencecondition is set as the Euclidean length of the mean shift, vector issmaller than a preset threshold.

Actually, in the mean shift based algorithm as used in these teachings,the mean shift update thread is performed multiple times by taking eachpixel in the image plane as the starting point and replace the currentpixel with the converged local maxima point. All the pixels leading tothe same local maxima will be set as the same label in the label image.After this, the very first mean shift segmentation (strictly speaking,it is the mean shift smooth filtering) result is obtained while it isalmost definitely over-segmented. Therefore, the over-segmented imagehas to be merged based on some rules. In the fusion step, extensive usewas made of region adjacency graphs (RAG).

The method flowchart is shown as in FIG. 2a . The reason for using theLUV color space is because that perceived color differences shouldcorrespond to Euclidean distances in the color space chosen to representthe features (pixels). The LUV and LAB color space were especiallydesigned to best approximate uniformly color space. To detect all thesignificant modes, the following basic mean shift filtering processshould be run for multiple times (evolving in principle in parallel)with different starting points that cover the entire feature space. Inthese teachings, all the pixels in the image domain are used as thestarting points.

The first step in the mean shift based feature space with the underlyingdensity ƒ(x) is to find the modes of this density. The modes are locatedamong the zeros of the gradient ∇ƒ(x)=0, and the mean shift procedure isan elegant way to locate these zeros without estimating the densitycompletely. The mean shift vector m(x) computed in Eq. 3.6 with kernel gis proportional to the normalized density gradient estimate obtainedwith kernel k. In Eq. 3.6, “n” represents the number of neighbor pixelsx_(i) involved in the kernel density estimation see, for example, C. M.Christoudias, B. Georgescu, P. Meer, Synergism in Low Level Vision, IEEEProc. of 16^(th) Inter. Conf. on Pattern Recognition, 2002, Vol. 4: pp.150-155, which is incorporated by reference herein is entirety and forall purposes). All the neighbor pixels located within the Euclideandistance h from the current pixel will be chosen. The mean shift vectorthus always points toward the direction of the maximum increase in thedensity. In this case, the y_(k) is iteratively updated according to Eq.3.7 until the convergence will lead to the local maxima for the currentpoint in the probability density function (PDF). The convergence isdefined as when the difference between y_(k) and y_(k+1) is smaller thana specified threshold value.

In Eq. 3.6, “i” represents the i^(th) gradient descent path. After allthe local maxima have been detected from different starting points, allthe points on the path leading to each maxima will be claimed to belongto the basin marked by the current maxima. Then the basins with the sizesmaller than the pre-stetting threshold value will be merged to thenearest basin whose size is bigger than a preset threshold. In bothequations, the pixel is described by a 5 dimension vector concatenatedin the joint spatial-range domain including 3 elements for the LUV colordomain and 2 elements for the spatial domain. As stated hereinabove, thekernel function k is chosen as the Epanechinikov kernel. In theseteachings, the combined kernel function shown in Eq. 3.8 is used.

$\begin{matrix}{{K_{{hs},{hr}}(x)} = {\frac{C}{{hs}^{2}{hr}^{3}}{k\left( {{\frac{x^{s}}{h^{s}}}}^{2} \right)}{k\left( {{\frac{x^{r}}{h^{r}}}}^{2} \right)}}} & {{Eq}.\mspace{14mu} 3.8}\end{matrix}$where hs and hr are different bandwidth values for spatial domain andrange domain, respectively.

After the initial mean shift filtering procedure, the over-segmentedimage are merged based on some rules. In the fusion step, extensive usewas made of region adjacency graphs (RAG). The initial RAG was builtfrom the initial over segmented image, the modes being the vertices ofthe graph and the edges were defined based on 4-connectivity on thelattice. The fusion was performed as a transitive closure operation onthe graph, under the condition that the color difference between twoadjacent nodes should not exceed hr/2.

Mean Shift Based Segmentation

The wound boundary determination approach using mean shift basedsegmentation is theoretically full-automatic and does not depend on anya priori manual input, which makes it computationally very economic andflexible for implementation in any hardware platforms. In FIG. 2b , theworkflow of this approach is provided.

The mean shift based algorithm is first applied to segment the originalwound image into a number of homogeneous regions. The mean shiftalgorithm is chosen over other segmentation methods, such as level setand graph cut based algorithms for several reasons. First, the meanshift algorithm takes into consideration the spatial continuity insidethe image by expanding the original 3D color range space to 5D space,including two spatial components, since direct classification on thepixels proved to be inefficient. Second, a number of accelerationalgorithms are available. Third, for both mean shift filtering andregion merge methods, the quality of the segmentation is easilycontrolled by the spatial and color range resolution parameters. Hence,the segmentation algorithm can be adjustable to different degrees ofskin color smoothness by changing the resolution parameters. Finally,the mean shift filtering algorithm is suitable for parallelimplementation since the basic processing unit is the pixel. In thiscase, the high computational efficiency of GPUs can be exploited, whichcan further improve the efficiency and achieve the real time woundassessment even on the smartphone-alone system.

After applying the mean shift algorithm, the image is usuallyover-segmented, which means that there are more regions in thesegmentation result than necessary for wound boundary determination. Tosolve this problem, the over-segmented image is merged into a smallernumber of regions which are more object-representative based on somerules. In the fusion step, extensive use was made of region adjacencygraphs (RAG). The initial RAG was built from the initial over-segmentedimage, the modes being the vertices of the graph and the edges definedbased on 4-connectivity on the lattice. The fusion was performed as atransitive closure operation on the graph, under the condition that thecolor difference between two adjacent nodes should not exceed h_(ƒ),which is regarded as the region fusion resolution. Based on experimentalresults, the over-segmentation problem is found to be effectively solvedby region fusion procedure.

Foot Outline Detection and Task Categorization

The wound boundary determination method is based on two assumptions.First, the foot image contains little information not related to thechronic wound. In reality, it is not a critical problem as it is assumedthat the patients and/or caregivers will observe the foot image with thewound on the smartphone screen before the image is captured to ensurethat the wound is clearly visible. Second, it is assumed that thehealthy skin on the sole of the foot is a nearly uniform color feature.

The largest connected component detection is first performed on thesegmented image, using the fast largest connected component detectionmethod including two passes. In the processing step Foot ColorThresholding, the color feature, extracted in the mean shiftsegmentation algorithm of this component, is compared with an empiricalskin color feature by calculating the Euclidean distance between thecolor vector for the current component and the standard skin colorvector from the Macbeth color checker. If the distance is smaller than apre-specified and empirically determined threshold value, the foot areais considered as having been located. Otherwise, the largest componentdetection algorithm is iteratively repeated on the remaining part of theimage while excluding the previously detected components until the colorthreshold condition is satisfied. After the foot area is located, abinary image is generated with pixels that are part of the foot labeled“1” (white) and the rest part of the image labeled “0” (black).

Then, the wound boundary determination tasks have to be classified intotwo categories: 1) the wound is fully enclosed within the foot outline;2) the wound is located at (or very near to) the boundary of the footoutline. The initial idea was to use the foot boundary smoothness todistinguish between these two situations. However, the problem is that agold standard for the ordinary smooth foot curve may be needed, i.e. theboundary of the healthy foot, and quantitatively compare the actuallydetected foot outline to it in some way. The search for such a groundtruth healthy foot curve is never an easy task. Moreover, it has to beensured that the patient's entire foot is imaged completely, which is adifficult-to-meet expectation for a self-management wound analysissystem considering the possible low mobility and the lack of experienceof handheld communication device use for most type 2 diabetic patients.Therefore, the following method is used to realize the taskclassification.

At first, one of the image morphology operations called a closingoperation (with a 9×9 circle structure element) is applied to remove allthe holes in the foot region (white part in the binary image) and smooththe external foot boundary, which will help us to eliminate the possibleinterference for accurate wound boundary determination. Secondly, thecombined region and boundary algorithm is applied to trace the externalfoot boundary along the edge of the white part in the foot binary image,as well as all the internal boundaries if there are any. For all theinternal boundaries in a foot region, only the ones with the perimeterlarger than a preset threshold (in one embodiment, it is set as 50 pixellengths) are kept. This simple threshold method may not be a perfectalgorithm but it works for most of the wound images in many instances.In other words, if there is at least one internal boundary exceeding thepreset threshold within the foot region, it is regarded as the woundboundary and returns it as the final boundary determination result. Onthe other hand, if there are not any internal boundaries whose lengthare beyond the threshold, other boundary determination methods may beneeded. Note that here it is assumed there is at least one wound area onthe photographed foot.

Boundary Determination for the Wounds Near the Edge of the Foot Outline

After a careful study and observation, a wound boundary determinationmethod, as shown in the right column in FIG. 2b , is used that isapplicable for a wound located at or near to the foot outline. As statedherein above, the external boundary of the non-enclosed foot outline isalready available.

As illustrated in the block diagram in FIG. 2b , the input to thismethod is the external boundary. Instead of keeping all the points onthe foot boundary, the edge points (all the points on the externalboundary of the foot region) are down-sampled by applying the HarrisCorner Detection method to a number of corner points (as shown in part(a) of FIG. 3). The corner points, also called junctions of edges, areprominent structural elements in an image and are therefore useful in awide variety of computer vision application. It is claimed that cornerpoints are more robust features for geometric shape detection than theregular edge points. In one instance, the perimeter of the foot outlineusually is made up of over 2000 pixels. After down sampling by cornerdetection, the number of corner points is approximately around 60-80(also in terms of pixels). This will greatly improve the time efficiencyof the algorithm. Besides, this down-sampling procedure will alsobenefit us when detecting the turning points (this will be discussed indetail in herein below).

The third to the eighth blocks in the right column in FIG. 2b shows themain idea, which is to detect the three turning points on the footboundary. These turning points will be used to determine the woundsection on the foot outline (shown in the part (c) of FIG. 3 by markingthe three turning points with small black crosses). The turning pointscan be defined as the most abruptly changing points along the footboundary. After the three points are determined, one can move along fromthe global maximum point (which is the most concave point in the middle)to the two local minimum points along the foot boundary in two oppositedirections. Then the two local minimum points are connected by an arcwhich is the substantially optimal approximation to the non-closed partof the wound boundary. In one instance, an arc is drawn from eitherlocal minimum to the other one with a radius equal to half of thediagonal length for the wound image. The trace generated by the abovetwo steps is the target wound boundary (as shown by a red curve in FIG.3d ).

For detecting the turning points, a maximum-minimum searching approachis used to detect the turning points. Herein below, a detaileddescription of this approach is provided.

First, all the corner points are sorted into a list based on theirposition on the foot boundary (from the top-right to top-left, in aclock-wise direction), then locate the two special extreme corner pointson the foot boundary: the leftmost and the rightmost (as indicated inpart (a) of FIG. 3). Then, the corner points are divided into twogroups: the corners points which located between the two extreme pointsand the corner points located outside this range. Note that thiscategorization is based on the clockwise sorted corner points list. InFIG. 4a , the first group of corner points is marked by U and the secondgroup is marked by L. For the first group the vertical distance of eachcorner point to the top side of the Smallest Fitting Rectangle (SFR) ofthe foot region is calculated. The smallest fitting rectangle issupposed to be tangent to the foot area at four boundary points: thetop-most, bottom-most, left-most and right most point of the foot area,as shown by the frame in FIG. 4 b.

Similarly, the vertical distance of each corner point in the secondgroup to the bottom side of the SFR (as shown in FIG. 4b ) iscalculated. Afterwards, the turning points are located by seeking forthe corner point with global maximum vertical distance to thecorresponding rectangle side (top or bottom based on its group number:first or second) and also two corner points on each side of the maximumpoint with the smallest local minimum vertical distance (as shown FIG.3c ). The only concern is the search for target turning points may beaccidentally stopped by interfering local extrema, which is a commonproblem of most local search algorithms. As stated above, a certainnumber of corner points are kept on the foot outline. Based onexperimental results, it is found that this down-sampling procedure caneliminate most of the interfering local extrema which may impede thesearch for the optimal turning points.

The above disclosed method mainly classifies the wound locations intothree categories: 1) wound in the middle of the foot, 2) wound at theedge of the foot without toe-amputation and 3) wound at the edge of thefoot with toe-amputation. For the first category, the wound is supposedto be surrounded by healthy skin and can be easily detected by tracingthe internal boundary within the foot outline. For the second and thirdcategories, the three turning points detection method is applied tolocate the wound boundary which is assumed to be the most irregularlychanged section on the foot outline. In practice, the method dealingwith these two situations (with or without toe-amputation) is slightlydifferent. Hence, the toe-amputation information may need to be given asan input to the method and obtained as part of the patient's medicalinformation.

In another instance, in the method of these teachings, analyzing theimage includes using a trained classifier and, in the system of theseteachings, the image analysis component is configured to use a trainedclassifier.

A machine learning based solutions has been developed in which the woundboundary determination is an object recognition task since it is claimedthat the machine learning (ML) is currently the only known way todevelop computer vision systems that are robust and easily reusable indifferent environments. Herein below, the term “wound recognition” isused as the equivalent expression of “wound boundary determination”,since both have the same goal.

In object recognition field, three major tasks needed to be solved toachieve the best recognition performance: 1) find the bestrepresentation to distinguish the object and background, 2) find themost efficient object search method and 3) design the most effectivemachine learning based classifier to determine whether a representationbelongs to the object category or not.

For the Chronic wound recognition method and components of theseteachings, a hybrid of the global window and local patch basedrepresentation which modifies the general form of the global texturedescriptors to be extracted within only local sub-windows or patches isused. A popular approach of this hybrid type is called Bags of VisualWords (BoW) which uses a visual vocabulary to compactly summarize thelocal patch descriptors within a region using a simple 1D histogram(see, for example; (i) Fei-Fei Li; Perona, P., “A Bayesian HierarchicalModel for Learning Natural Scene Categories”. 2005 IEEE Computer SocietyConference on Computer Vision and Pattern Recognition (CVPR'05) p. 524and (ii) Rob Fergus, Classical Methods for Object Recognition, slidespresented at ICCV 2009 course, both of which are Incorporated byreference herein in their entirety and for all purposes).

This representation is completely orderless, which means that greaterflexibility is allowed (for better or worse) with respect to viewpointand pose changes. At the same time, the invariance properties of theindividual local descriptors make them a powerful tool to tolerate thevariation of the viewpoint or pose while giving informative localappearance cues. The regularity or rigidity of an object category'sappearance pattern in 2D determines which style is better suited. Forexample, the class of frontal face is quite regular and similaritystructured across instances, and thus is more suitable for the 2Dlayout-preserving descriptors; in contrast, wounds represent a varietyof shapes of which most are irregular. This property makes it suited toa more flexible summary of the texture and key features. What isparticularly convenient about the bag-of-words (BOW) representation isthat it translates a (usually very large) set of high-dimensional localdescriptors into a single sparse vector of fixed dimensionality acrossall images. This in turn allows one to use many machine learningalgorithms that by default assume that the input space isvectorial—whether for supervised classification, feature selection, orunsupervised image clustering.

The object localization techniques can mainly fall into one of twocategories: 1) the “top-down” technique, which tries to fit a coarseglobal object model to each possible location on the image grid or 2)the “bottom-up” technique, which tries to produce a pixel levelsegmentation of the input image and are built from the bottom up onlearned local representation and can be seen as an evolution of texturedetectors. The sliding window technique is a typical example for thefirst category. Due to its algorithmic nature, the sliding window searchapproach suffers from several limitations including the highcomputational cost, little room for error and inflexible for accuratewound boundary determination. Hence, in some related references, it isclaimed the bottom-up technique is more suitable for object classsegmentation task (similar as the wound boundary determination task).

The supervised learning methods have been most widely used. Thisapproach will try to inferring a model from labeled training data. Inrelated references, the comparison of several most popular supervisedlearning methods is provided. Generally speaking, support vector machine(SVMs) tends to perform much better when dealing with multi-dimensionsand continuous features. (See, for example, Using Support Vector Machineas a Binary Classifier, International Conference on Computer Systems andTechnologies—CompSysTech' 2005 which is incorporated by reference hereinin its entirety and for all purposes). For SVMs, given a set of trainingexamples, each marked as belonging to one of two categories, an SVMtraining algorithm builds a model that assigns new examples into onecategory or the other. An SVM model is a representation of the examplesas points in space, mapped so that the examples of the separatecategories are divided by a clear gap that is as wide as possible. Newexamples are then mapped into that same space and predicted to belong toa category based on which side of the gap they fall on. In addition toperforming linear classification, SVMs can efficiently perform anon-linear classification using what is called the kernel trick, whichimplicitly mapping their linear inputs into high dimensional featurespaces.

In one embodiment of these teachings, a two-stage recognition schemebased on some object recognition approach already being successfullyapplied in the pedestrian recognition task is used. The workflow of thiswound recognition system is shown in FIG. 5.

In the training process, there are two stages. For both stages, the SVMbased binary classifier training method is used. In the first stage, thesuper-pixel segmentation is performed by either the mean shift or SLIC(simply linear iterative clustering) algorithm to group pixels intoperceptually meaningful atomic regions which can be used to replace therigid structure of the pixel grid. Then, the vector representation isbuilt up for each super-pixel by using the bag of words (BCW) histogrambased on local DSIFT (dense SIFT) or SURF feature descriptor within thecurrent super-pixel.

To generate this representation, the extracted descriptors are thenquantized using a K-means dictionary and aggregated into one normalizedhistogram h_(i)∈R₊ ^(K) for each super-pixel s_(i) in the image, where Kis the number of words predefined in the codebook (the set of clustersresulted from the K-means algorithm). In order to train a classifier,each super-pixel s_(i) is assigned the most frequent class label itcontains (in this case, some manually labeled ground truth images whichhave pixel-level granularity are needed). Then a SVM with an RBF kernelis trained on the labeled histograms for either category: wound andnon-wound. This yields discriminant functions is proposed in relativereferences and shown as below.

$\begin{matrix}{{C(h)} = {\sum\limits_{j = 1}^{L}\;{c_{i}{\exp\left( {{- \gamma}\;{d^{2}\left( {h,h_{i}} \right)}} \right)}}}} & (1)\end{matrix}$where c_(i)∈R are coefficients and h_(i) representative histograms(support vectors) selected by SVM training, γ∈R⁺ is a parameter selectedby cross-validation, and d²(h, h_(i)) is the vector distance between thecurrent histogram h and each support vector.

This classifier which results from this is very specific. It findssuper-pixels which resemble super-pixels that were seen in the trainingdata without considering the surrounding region. However, a drawback oftraining a classifier for each super-pixel is that the histogramsassociated with each super-pixel are very sparse, often containing onlya handful of nonzero-elements. This is due to the nature of thesuper-pixels: by definition they cover areas that are roughly similar incolor and texture. Since the features are fixed-scale and extracteddensely, the super-pixels sometimes contain tens or even hundreds ofdescriptors that quantize to the same visual word.

To overcome the problems caused by the lack of consideration of thesurrounding region of each super-pixel and sparse histogramrepresentation, the histograms are applied based on super-pixelneighborhoods. Let G(S, E) be the adjacency graph of super-pixels s_(i)in an image, and h_(i) ⁰ be the non-normalized histogram associated withthis region. E is the set of edges formed between pairs of adjacentsuper-pixels (s_(i), s_(j)) in the image and D(s_(i), s_(j)) is thelength of shortest path between two super-pixels. Then, h_(i) ^(N) isthe histogram obtained by merging the histograms of the super-pixels_(i) and neighbors who are less than N nodes away in the graph:

$\begin{matrix}{h_{i}^{N} = {\sum\limits_{S_{j}❘{{D{({s_{i},s_{j}})}} \leq N}}\; h_{j}^{0}}} & (2)\end{matrix}$The training framework is unchanged, except that super-pixels aredescribed by the normalized histograms h_(i) ^(N) in place of h_(i).

Finally, these 1D merged histogram representations are taken as theinput for the binary SVM training module. After the binary classifier istrained, it is applied to classify all super-pixels from all trainingimages. Then, all the super-pixels labeled as wound are gathered by thefirst stage classifier and an approximately equal number of non-woundsuper-pixels as the training data set for the next stage of machinelearning. For each instance in this set, the dominant color descriptor(DCD) is extracted and train the second stage classifier (whichinherently shares the same working scheme with the first stage SVM basedclassifier) based on these descriptors.

In order to compute this descriptor, the colors present in a givenregion are first clustered. This results in a small number of colors andthe percentages of these colors are calculated. As an option, thevariances of the colors assigned to a given dominant color are alsocomputed. The percentages of the colors present in the region should addup to 1. A spatial coherency value is also computed that differentiatesbetween large color blobs versus colors that are spread all over theimage. The descriptor is thus defined as following:F={(c _(i) ,p _(i) ,v _(i)),s},(i=1,2, . . . ,N)  (3)where c_(i) is the i^(th) dominant color and p_(i) is its percentagevalue and v_(i) is its color variance. N represents the number ofdominant color clusters. The spatial coherency s is a single number thatrepresents the overall spatial homogeneity of the dominant colors in theimage. In one instance, the DCD can be easily determined from the earlymean shift based super-pixel segmentation results. The reason for thesecond stage classification is to utilize the color features to furtherimproving the differentiation between skin and wound tissues near thewound boundary.

In the testing process, for an input testing image, same super-pixelsegmentation and BoW representation generation will be performed. Then,the first stage binary classifier is applied to identify all “candidatewound” super-pixels. Next, the DCD descriptor is generated for each“candidate wound” super-pixel and input to the second stage binaryclassifier. Next, a conditional random field (CRF) technique basedrefinement method is operated to recover more precise boundaries whilestill maintaining the benefits of histogram merge over the super-pixelneighborhood. Finally, a closing operation, one of the morphologymethods, can be performed to eliminate small holes in the detected woundarea and further to smooth the wound boundary.

To train the classifier and also evaluate the wound recognitionperformance of the method of these teachings, the help of experiencedwound clinicians is needed to generate the ground truth wound labels. Inone instance, 48 wound images collected from UMass Wound Clinic from 12patients over 12 months are used. For each image, three clinicians wereasked to delineate the wound boundary independently with Photoshopsoftware and a set of electronic drawing pen and panel. Afterwards, themajority vote scheme is used (for each pixel, if 2 or 3 clinicians labelit as “wound”, then it will be determined as “wound” pixel. Otherwise,it will be determined as “non-wound” pixel). An example of the groundtruth generation is illustrated in FIG. 6.

The samples of the wound recognition results on the images of realpatients are shown in FIG. 7. It can be seen that this solution providepromising wound boundary determination.

In order to better assess the SVM based wound recognition method, thefollowing testing and evaluation approach is used. First, theleave-one-out cross validation method is adopted to evaluate the modelperformance on the entire dataset. Specifically, one image is choseneach time from the sample image set as the testing sample and the restis taken as the training samples used for SVM based model training.Hence, this experiment has to performed for a number of times equal tothe size of the entire sample image set (48 times for all 48 woundimages) in order to test on the entire image dataset and keep thespecified testing image different from all images in the trainingdataset.

Second, since the wound recognition is a skewed distributed binary classproblem, which contains a large number of non-wound super-pixels and arelatively small number of wound super-pixels for each wound image, theaccuracy rate cannot be used to evaluate the performance. Instead, theidea of true positive (tp), false positive (fp), false negative (fn) andtrue negative (tn) respectively defined as in FIG. 8 is used. The matrixshown in this figure is also called the confusion matrix.

Substantial research has been performed on developing a convincingevaluation score based on these four values. The Matthews CorrelationCoefficient (MCC) is used in machine learning as a measure of thequality of binary classification. Especially, it takes into account trueand false positives and negatives and is generally regarded as abalanced measure which can be used even if the classes are of verydifferent sizes.

The MCC is in essence a correlation coefficient between the observed andpredicted binary classification; it returns a value between −1 and +1. Acoefficient of +1 represents a perfect prediction, 0 no better thanrandom prediction and −1 indicates total disagreement between predictionand observation. It is defined directly on the confusion matrix asbelow.

$\begin{matrix}{{MCC} = \frac{{{tp} \times {tn}} - {{fp} \times {fn}}}{\sqrt{\left( {{tp} + {fp}} \right)\left( {{tp} + {fn}} \right)\left( {{tn} + {fp}} \right)\left( {{tn} + {fn}} \right)}}} & (4)\end{matrix}$The experimental results shows that the average MCC value of 48 testimages using leave-one-out evaluation method is 0.7, which is 0.1 higherthan the commonly regarded standard value of promising objectrecognition.

In one instance, in the method of these teachings, performing the colorsegmentation comprises performing a K-mean color clustering algorithm;and evaluating the wound area comprises using a red-yellow-blackevaluation model for evaluation of the color segmentation and, in thesystem of these teachings, the image segmentation component isconfigured to perform a K-mean color clustering algorithm; and uses ared-yellow-black evaluation model for evaluation of the colorsegmentation.

Red-Yellow-Black (RYB) Model

After the accurate wound boundary is acquired, the wound area isanalyzed within the boundary using some wound description model. Manymethods for assessing and classifying open wounds require advancedclinical expertise and experience. Specialized criteria have beendeveloped for diabetic foot ulcers. In order to facilitate the woundmanagement performed by patients themselves at home, there is need for asimple classification system that can be universally applied. The RYBwound classification model which was first proposed in the October 1988by J. Z. Cuzzell and C. Blanco provide us a consistent, simple model toevaluate the wound (D. Kransner, Wound Care How to Use theRed-Yellow-Black System, the American Journal of Nursing, Vol. 95 (5),1995, pp. 44-47 which is incorporated by reference herein in itsentirety and for all purposes).

The RYB system classifies the wound as red, yellow, black or mixedtissues which represent the different phases of the tissue on thecontinuum of the wound healing process, respectively. In detail, redtissues are viewed as the inflammatory (reaction) phase, proliferation(regeneration), or maturation (remodeling) phase. On the other hand,yellow tissues stand for the infected or contain slough that aren'tready to heal. At last, black tissues indicate necrotic tissue state,which is not ready to heal either.

Based on the RYB wound evaluation model, the task for wound analysis isequal to clustering all the pixels within the wound boundary intocertain color categories. Therefore, all classical clustering method canbe applied to solve this task.

K-Mean Algorithm

In data mining, k-means clustering is a method of cluster analysis (see,or example, K means and Hierarchical Clustering, tutorial slides byAndrew W. Moore, 2001, which is incorporated by reference herein in itsentirety and for all purposes), which aims to partition n observationsinto k clusters in which each observation belongs to the cluster withthe nearest mean. In one instance, all the pixels within the woundboundary can be viewed as observations. The three colors referred in RYBmodel are regarded as clusters. The algorithm is graphically illustratedin FIG. 9 a.

There are several things needed to be further specified.

1) The color difference between the cluster center and the target pixel(expressed as Eu in the flowchart in part a) in FIG. 9c is calculated bythe standard Euclidean color difference in CIE Lab model.

2) Strictly speaking, K-mean algorithm is a NP-hard problem, which isunable to converge to a solution within limited time when the image sizeis large enough. However, the iteration can be terminated when theaverage mean variance of each cluster is smaller than a pre-specifiedthreshold. This heuristic method is expressed as the decision block inpart a) of FIG. 9c . In part a) of FIG. 9c , the initial centers arechosen randomly. However, in practice, the initial centers may bespecified according to some empirical values such as the Macbeth ColorChecker. By this operation, the converging speed will be increased thusmaking the color clustering process more efficient.

3) As shown in FIG. 9a , the number of cluster is preset to 3. However,the number could be smaller or larger than 3. Some post-processing hasto be performed to the resulting clusters. In the present teachings,only the situation that the number of clusters is equal to 3 at most isconsidered. Therefore, some two or three clusters may have to becombined if they are close to each other enough, which can be equallyviewed as the mean value difference of the two or three clusters issmaller than a preset threshold.

In another embodiment, in the method of these teachings, performing thecolor segmentation includes using a K-mean color clustering algorithm onresults of images used on training a classifier and, in the system ofthese teachings, the image segmentation component is configured to use aK-mean color clustering algorithm on results of images used on traininga classifier.

A method similar to Bag-of-Words is used in another embodiment of colorsegmentation of these teachings. The flow chart of this embodiment isshown in FIG. 9b . There are two major tasks in this algorithm. In thetraining process, as stated herein above, three experienced woundclinicians were asked to label the wound area using a set of electronicdrawing panel, pen and also Photoshop software on the laptop. First, thecolor vectors are gathered, in CIE Lab color space, for all labeledwound pixels in the sample foot ulcer images. Then, the K-meanclustering algorithm are performed for all color vectors. Instead ofusing the pre-set standard color center vector, a number of color vectorvalues are randomly selected as the initial centers. It turns out thatthe setting of the initial centers has no impact on the final colorclustering results. The only parameter, which needs to be preset is thecluster number. In one instance, it is set to 10, a relatively largecluster number considering the wound tissues usually only contain 3-5obviously distinct color types, since this will provide us a more fineclassification and also is not too time demanding, which means eachpixel is assigned to a cluster center reasonably resembling its owncolor. After the initial clustering, all the cluster centers areanalyzed and several centers with small Euclidean distance in the colorspace are merged into one. This operation can reduce the final clusternumber and form a more representative wound evaluation model. From theten color cluster centers resulted from the K-mean algorithm, only anumber of colors centers (1, 2, 3, 4 in one instance) are quite distinctfrom each other. Those color centers can be regarded as clusters foryellow, white, black and red, respectively. Some of the other colorcenters (in one instance from 5-9) can all be classified as one color(red, in one instance) but with different saturation and hues or mergedinto one of the distinct color centers. The original RYB model isextended to include another cluster center representing the wound tissuein color of white. However, based on clinicians' opinion, the whitewound tissue, which is more like the calluses, should not contribute tothe final wound healing assessment. Hence, the white tissue isconsidered as part of the wound area but is not considered whenperforming the color based wound analysis.

In the segmentation process, the original set of clusters (in number of10) is used and the assignment is made to each wound pixel in thedetermined wound area. After that, the pixels assigned to cluster number5-9 are merged into cluster 4, and the pixels assigned to cluster number10 are merged into cluster 3, since the Euclidean distance in CIE Labcolor space is small enough. The color segmentation results on 5 samplewound images are shown in FIGS. 10a-10i , the original images aredisplayed in 10 a-10 c, and the color segmentation results using theK-mean algorithm alone shown in FIGS. 10d-10f . The results from theabove described color segmentation algorithm in this report can be seenin FIGS. 10g-10i . After comparison, algorithm combining k-means withthe figures selected by the clinicians results in an improvement.

In one instance, in the method of these teachings, evaluating the woundarea includes determining a healing score as a method for quantifying ahealing status of the wound area and, in the system of these teachings,the wound evaluation component is configured to determine a healingscore as a method for quantifying a healing status of the wound area.

Healing Score

One goal of these teachings is to provide more meaningful wound analysisresults to the users, including both the clinicians and diabeticpatients. For clinicians, the wound area size and different color tissuecomposition may be sufficient. They can make their diagnosis based onthese raw data. However, for ordinary patients assumed to be without anyclinical knowledge about wounds, only providing them some absolutenumbers does not give them with much help in understanding their actualwound status. Hence, there is a need to translate the raw data into ameaningful numerical value, like a score in the range of 0-100, wherelarger simply means better. In this report, a numerical wound evaluationvalue called healing score is used. The basis for calculating thehealing score are four indicators: wound area size, red tissue size,yellow tissue size, and black tissue size. As introduced in relatedreferences, the red means granulation, which is probably a positive signfor healing. On the other hand, yellow might represent tissues withinfection and black stands for necrotic tissues. And these are negativesigns for bad healing status. Besides, the shrinking of the entire woundarea certainly is a strong positive evidence of good healing status.Note that since there is no official described clinical correspondencefor the white tissue, only the red, yellow and black tissues areconsidered for the calculation of the healing score and will merge thewhite cluster to the closet one of the three types.

Considering all of the factors above, a healing score calculationformula is provided herein below. The Healing Score formulation hasthree components:

1) A Healing Score based on wound area, which will have an initial scoreof 50, given that the wound area can become larger or smaller

2) A Healing Score based on the color with the wound boundary. Here, theinitial score is not fixed, but will be bounded by the range 0-100, suchthat all red will produce a Healing Score of 100, all black will producea Healing Score 0, and some combination of red, white, yellow and blackwill generate a Healing Score 0<score<100.

3) A composite Healing Score, which will be a weighted average of theHealing Score based on wound area and the Healing Score based on thecolor with the wound boundary. The weight may be constant or may beinfluenced by the size of the wound area.

Healing Score Based on Wound Area

As stated, the initial value is defined to be 50. Let α_(n) be the woundarea in week n and S_(n) ^(A) be the wound area score in week n. a₀ isthe initial wound area size acquired when the patient use the system forthe first time. Thus, S_(n) ^(A)=ƒ(a₀,a_(n)) and S₀ ^(A)=ƒ(a₀,a₀)=50,where f is supposed to be function taking a_(n) and a₀ as itsparameters.

$\begin{matrix}{{{S_{n}^{A} = {\left( {1 - \frac{a_{n} - a_{0}}{a_{0}}} \right)50}},{a_{n} \leq {2a_{0}}}}{S_{n}^{A} = 0},{a_{n} > {2a_{0}}}} & (5)\end{matrix}$As a_(n) varies from 0 to 2a₀, S_(n) ^(A) decreases linearly from 100 to0. For values of an a_(n)>2a₀, S_(n) ^(A)=0. This should be reasonableassumption that once the wound become twice as large as the initial sizethere is no sign of healing at all. The wound area healing score is arelative numerical value which takes the initial wound area size as thereference.

Healing Score Based on Color with the Wound Boundary

Let S_(n) ^(T) be the Healing Score based on the color with the woundboundary in week n. Similar to a_(n), the ratio of red area, yellow areaand black area are defined, within the wound boundary, as r_(n), y_(n)and b_(n), respectively, and where subscript ‘n’ refers to week n.Clearly, r_(n)+y_(n)+b_(n)=1 in general, and specifically r₀+y₀+b₀=1.

Based on wound evaluation theory, S_(n) ^(T) must be formulated so thatS_(n) ^(T)=100 for r_(n)=1; y_(n)=b_(n)=0, and S_(n) ^(T)=0 for b_(n)=1;r_(n)=y_(n)=0. The following formulation for S_(n) ^(T) is proposed:

$\begin{matrix}{S_{n}^{T} = {\frac{1 + r_{n} - {0.5y_{n}} - b_{n}}{2}100}} & (6)\end{matrix}$It is easily verified that S_(n) ^(T)(r_(n)=1; y_(n)=b_(n)=0)=100 andthat S_(n) ^(T)(b_(n)=1; y_(n)=r_(n)=0)=0. Consider also the case wherer_(n)=y_(n)=b_(n)=0333, giving S_(n) ^(T)=41.7.

Composite Healing Score

Let S_(n) be the overall, or composite, Healing Score:S _(n) =w _(A) S _(n) ^(A) +w _(T) S _(n) ^(T)  (7)where w_(A) and w_(T) are weights, such that w_(A)+w_(T)=1. This allowsus to formulate S_(n) asS _(n) =w _(A) S _(n) ^(A)+(1−w _(A))S _(n) ^(T)  (8)A simple (and acceptable) solution is to set w_(A)=0.5. w_(A) does nothave be a constant; instead, w_(A) should have a greater influence whenthe wound is close to being healed and hence the area is small.Specifically, in one instance, w_(A) increases linear from w_(A)=0.5 tow_(A)=1.0, as S_(n) ^(A) increases linearly from 0 to 100. In otherwords,

$\begin{matrix}{{{w_{A} = {0.5 + {\frac{0.5}{100}S_{n}^{A}}}},{giving}}{S_{n} = {{\left\lbrack {0.5 + {0.05S_{n}^{A}}} \right\rbrack S_{n}^{A}} + {\left\lbrack {0.5 - {0.005S_{n}^{A}}} \right\rbrack S_{n}^{T}}}}} & (9)\end{matrix}$

An example of applying the proposed healing score to evaluate the woundstatus is based on five images. The wound analysis data for these fiveimages are shown in Table 1. After calculation, the healing score forthese four wound images are 82.78, 87.86, 84.75, and 75.59 (the firstimage is viewed as the base reference and not scored). From Image 1 to2, the wound area is shrinking. From Image 2 to 3, only a small sizedecrease of the wound area is observed. Hence, there is also a tinyincrease of the healing score by 4.4 points. From part Image 3 to 4,more surgical sutures were exposed and more yellow tissues occurred. Onthe other hand, the size of the entire wound area didn't change toomuch. Corresponding to this trend, the healing score is nearly 3 pointslower than the previous time. Finally, from part Image 4 to 5, there areextra yellow tissues generated on the outer part of the wound and thered tissues are shrinking. On the other hand, the wound and black tissuearea are decreased in a tiny degree. Hence, the healing score decreasedby nearly 9 points.

TABLE 1 Wound assessment results (area unit: mm²) Image 1 Image 2 Image3 Image 4 Image 5 Healing score 82.78 87.86 84.75 75.59 Wound area1126.57 403.09 293.17 279.34 457.39 Red area 791.21 353.25 282.84 214.95106.46 Yellow area 246.42 39.22 10.33 43.41 324.31 Black area 88.9410.62 0 20.98 26.62

Image Capture Box

In one embodiment, the system of these teachings includes an imagingcomponent having a first front surface mirror and a second front surfacemirror, the second front surface mirror being disposed at a right angleto the first front surface mirror, the imaging component beingconfigured such that the body part is positioned above the first andsecond front surface mirrors and away from an axis bisecting the rightangle; and wherein the image acquisition device is positioned above thefirst and second front surface mirrors, away from the axis bisecting theright angle and on an opposite side of the axis bisecting the rightangle from the body part.

To ensure consistent image capture conditions and also to facilitate aconvenient image capture process for patients with type 2 diabetes, animage capture device was designed in the shape of a box. This device isreferred to as “the Image Capture Box”. The image capture box wasdesigned as a compact, rugged and inexpensive device that: (i) allowspatients to both view the sole of their foot on the screen of a devicehaving an image capture components (for example, a handheld portablecommunication device such as, but not limited to, a smartphone) and tocapture an image since the majority of patients' wounds occur on thesoles of their feet, (ii) allows patients to rest their feetcomfortably, without requiring angling of the foot or the image capturecomponent 135 (in one instance, a smartphone camera), as patients may beoverweight and have reduced mobility, and (iii) accommodates imageviewing and capture of left foot sole as well as right foot sole. Toachieve these objectives, two front surface mirrors 115, 125 are used,placed at an angle of 90° with respect to each other, and with thecommon line of contact tilted 45° with respect to horizontal. Aschematic drawing of basic optical principle for foot imaging is shownin FIG. 11. The optical path is represented by straight lines witharrows indicating the direction.

A SolidWorks™ 3D rendering of the image capture box is shown in FIGS.12a-12c . As seen in this figure, the entire box has a rectangulartrapezoid shape. Rectangular openings for placing the foot andsmartphone are cut into the slanted surface, shown in FIG. 12b , whichis at 45° with respect to horizontal. In this case, the patient can resthis/her foot comfortably and view his/her wound on the LCD display ofthe smartphone camera. When using the box, the patient needs to ensurethat the wound is completely located within the opening by simplyobserving the image displayed on the smartphone.

To avoid the ghost image effect associated with normal back surfacemirrors (reflective surface on the back side of the glass), frontsurface mirrors (reflective surface on the front side) are needed, asillustrated in FIG. 12a . The optical paths for both the front surfacemirror and the normal mirror are shown in FIG. 13.

In one embodiment, the image acquisition device, the image analysiscomponent, the image segmentation component and the wound evaluationcomponent of the system of these teachings are comprised in a handheldportable electronic device. In that embodiment, the handheld portableelectronic/communication device includes the image acquisition device,the image acquisition device being configured for capturing an image ofa body part including a wound area, one or more processors, and computerusable media having computer readable code embodied therein that, whenexecuted by the one or more processors, causes the one or moreprocessors to extract a boundary of the wound area, perform colorsegmentation within the boundary of the wound area, wherein the woundarea is divided into a plurality of segments, each segment beingassociated with a color indicating a healing condition of the segmentand evaluate the wound area.

Descriptions of exemplary implementations of a mobile-based system canbe found in, for example, U.S. Publication No. 2012/0190947 to Chon etal., which is incorporated herein by reference in its entirety for allpurposes.

FIG. 14 is a block diagram representation of one embodiment of thesystem of these teachings. Referring to FIG. 14, in the embodiment showntherein, a mobile communication system 280 includes a processor 250 andone or more memories 260. In the embodiment shown in FIG. 14, a camera265, where the camera as an objective lens 267, can also supply thephysiological indicators signal to the mobile communication device 280.The one or more memories 260 have computer usable code embodied thereinthat causes the processor 250 to that causes the processor to extract aboundary of the wound area, perform color segmentation within theboundary of the wound area, wherein the wound area is divided into aplurality of segments, each segment being associated with a colorindicating a healing condition of the segment and evaluate the woundarea. In one or more instances, the computer readable code causes theprocessor 250 to perform the implement the methods describedhereinabove.

The one or more memories 260 represent one embodiment of computer usablemedia having computer readable code embodied therein that causes aprocessor to implement the methods of these teachings. Embodiments ofthe method of these teachings are described hereinabove and the computerreadable code can cause a processor to implement those embodiments.

In the embodiment shown in FIG. 14, the mobile communication device 280also includes an antenna 265 that enables communications through one ormore of a variety of wireless protocols or over wireless networks.

In another embodiment, in the system of these teachings, the imageacquisition device is comprised in a handheld portable electronicdevice; and the image analysis component, and the wound evaluationcomponent are comprised in a computing component. The handheld portableelectronic device, such as that shown in FIG. 14, includes the imageacquisition device, the image acquisition device being configured forcapturing an image of a body part including a wound area, one or moreprocessors, and computer usable media having computer readable codeembodied therein that, when executed by the one or more processors,causes the one or more processors to transmit the image to the computingcomponent.

The computing component could have a structure such as that shown inFIG. 15. Referring to FIG. 15, in the structure shown there in, one ormore processors 155 are operatively connected to an input component 160,which could receive the images transmitted by the handheld portableelectronic/communication device, and to computer usable media 165 thathas computer readable code embodied therein, which, when executed by theone or more processors 155, causes the one or more processors 155 toperform the method of these teachings. The one or more processors 155,the input component 160 and the computer usable media 165 areoperatively connected by means of a connection component 170.

In one instance, the computer readable code embodied in the computerusable media 165 of the computing component causes the one or moreprocessors 155 to receive the image from the handheld portableelectronic device, extract a boundary of the wound area, perform colorsegmentation within the boundary of the wound area, wherein the woundarea is divided into a plurality of segments, each segment beingassociated with a color indicating a healing condition of the segmentand evaluate the wound area.

An exemplary embodiment of the system including a handheld portableelectronic/communication device and a computing component (also referredto as a collaborative or cooperative system, is shown in FIG. 16.Referring to FIG. 16, in the exemplary embodiment shown therein, thehandheld portable electronic/communication device is a smart phone andthe computing device is a laptop. The smartphone will play the role ofclient in the communication scheme. It will accomplish the followingtask sequentially: 1) take the picture of the wound and save on thespecified directory on the SD card; 2) make request to the server andsend the JPEG file of wound image to the laptop by function “post( )” inan http library and 3) receive the analyzed image file sent by laptopand display it on the screen. The laptop will be viewed as the serverparty, which will 1) listen to the request of the client and get theimage file sent by the client by “dopost( )” function, 2) accomplish thewound image process to he received wound image and 3) send the JPEG fileof the processed image back to the smartphone (client) as a response.

The following is a disclosure by way of example of a device configuredto execute functions (hereinafter referred to as computing device) whichmay be used with the presently disclosed subject matter. The descriptionof the various components of a computing device is not intended torepresent any particular architecture or manner of interconnecting thecomponents. Other systems that have fewer or more components may also beused with the disclosed subject matter. A communication device mayconstitute a form of a computing device and may at least include acomputing device. The computing device may include an inter-connect(e.g., bus and system core logic), which can interconnect suchcomponents of a computing device to a data processing device, such as aprocessor(s) or microprocessor(s), or other form of partly or completelyprogrammable or pre-programmed device, e.g., hard wired and orapplication specific integrated circuit (“ASIC”) customized logiccircuitry, such as a controller or microcontroller, a digital signalprocessor, or any other form of device that can fetch instructions,operate on pre-loaded/pre-programmed instructions, and/or followedinstructions found in hard-wired or customized circuitry to carry outlogic operations that, together, perform steps of and whole processesand functionalities as described in the present disclosure.

In this description, various functions, functionalities and/oroperations may be described as being performed by or caused by softwareprogram code to simplify description. However, those skilled in the artwill recognize what is meant by such expressions is that the functionsresult from execution of the program code/instructions by a computingdevice as described above, e.g., including a processor, such as amicroprocessor, microcontroller, logic circuit or the like.Alternatively, or in combination, the functions and operations can beimplemented using special purpose circuitry, with or without softwareinstructions, such as using Application-Specific Integrated Circuit(ASIC) or Field-Programmable Gate Array (FPGA), which may beprogrammable, partly programmable or hard wired. The applicationspecific integrated circuit (“ASIC”) logic may be such as gate arrays orstandard cells, or the like, implementing customized logic bymetallization(s) interconnects of the base gate array ASIC architectureor selecting and providing metallization(s) interconnects betweenstandard cell functional blocks included in a manufacturer's library offunctional blocks, etc. Embodiments can thus be implemented usinghardwired circuitry without program software code/instructions, or incombination with circuitry using programmed software code/instructions.

Thus, the techniques are limited neither to any specific combination ofhardware circuitry and software, nor to any particular tangible sourcefor the instructions executed by the data processor(s) within thecomputing device. While sonic embodiments can be implemented in fullyfunctioning computers and computer systems, various embodiments arecapable of being distributed as a computing device including, e.g., avariety of forms and capable of being applied regardless of theparticular type of machine or tangible computer-readable media used toactually effect the performance of the functions and operations and/orthe distribution of the performance of the functions, functionalitiesand/or operations.

The interconnect may connect the data processing device to define logiccircuitry including memory. The interconnect may be internal to the dataprocessing device, such as coupling a microprocessor to on-board cachememory or external (to the microprocessor) memory such as main memory,or a disk drive or external to the computing device, such as a remotememory, a disc farm or other mass storage device, etc. Commerciallyavailable microprocessors, one or more of which could be a computingdevice or part of a computing device, include a PA-RISC seriesmicroprocessor from Hewlett-Packard Company, an 80×86 or Pentium seriesmicroprocessor from Intel Corporation, a PowerPC microprocessor fromIBM, a Sparc microprocessor from Sun Microsystems, Inc, or a 68xxxseries microprocessor from Motorola Corporation as examples.

The inter-connect in addition to interconnecting such asmicroprocessor(s) and memory may also interconnect such elements to adisplay controller and display device, and/or to other peripheraldevices such as input/output (I/O) devices, e.g., through aninput/output controller(s). Typical I/O devices can include a mouse, akeyboard(s), a modem(s), a network interface(s), printers, scanners,video cameras and other devices which are well known in the art. Theinter-connect may include one or more buses connected to one anotherthrough various bridges, controllers and/or adapters. In one embodimentthe I/O controller includes a USB (Universal Serial Bus) adapter forcontrolling USB peripherals, and/or an IEEE-1394 bus adapter forcontrolling IEEE-1394 peripherals.

The memory may include any tangible computer-readable media, which mayinclude but are not limited to recordable and non-recordable type mediasuch as volatile and non-volatile memory devices, such as volatile RAM(Random Access Memory), typically implemented as dynamic RAM (DRAM)which requires power continually in order to refresh or maintain thedata in the memory, and non-volatile ROM (Read Only Memory), and othertypes of non-volatile memory, such as a hard drive, flash memory,detachable memory stick, etc. Non-volatile memory typically may includea magnetic hard drive, a magnetic optical drive, or an optical drive(e.g., a DVD RAM, a CD ROM, a DVD or a CD), or other type of memorysystem which maintains data even after power is removed from the system.

A server could be made up of one or more computing devices. Servers canbe utilized, e.g., in a network to host a network database, computenecessary variables and information from information in the database(s),store and recover information from the database(s), track informationand variables, provide interfaces for uploading and downloadinginformation and variables, and/or sort or otherwise manipulateinformation and data from the database(s). In one embodiment a servercan be used in conjunction with other computing devices positionedlocally or remotely to perform certain calculations and other functionsas may be mentioned in the present application.

At least some aspects of the disclosed subject matter can be embodied,at least in part, utilizing programmed software code/instructions. Thatis, the functions, functionalities and/or operations techniques may becarried out in a computing device or other data processing system inresponse to its processor, such as a microprocessor, executing sequencesof instructions contained in a memory, such as ROM, volatile RAM,non-volatile memory, cache or a remote storage device. In general, theroutines executed to implement the embodiments of the disclosed subjectmatter may be implemented as part of an operating system or a specificapplication, component, program, object, module or sequence ofinstructions usually referred to as “computer programs,” or “software.”The computer programs typically comprise instructions stored at varioustimes in various tangible memory and storage devices in a computingdevice, such as in cache memory, main memory, internal or external diskdrives, and other remote storage devices, such as a disc farm, and whenread and executed by a processor(s) in the computing device, cause thecomputing device to perform a method(s), e.g., process and operationsteps to execute an element(s) as part of some aspect(s) of themethod(s) of the disclosed subject matter.

A tangible machine readable medium can be used to store software anddata that, when executed by a computing device, causes the computingdevice to perform a method(s) as may be recited in one or moreaccompanying claims defining the disclosed subject matter. The tangiblemachine readable medium may include storage of the executable softwareprogram code/instructions and data in various tangible locations,including for example ROM, volatile RAM, non-volatile memory and/orcache. Portions of this program software code/instructions and/or datamay be stored in any one of these storage devices. Further, the programsoftware code/instructions can be obtained from remote storage,including, e.g., through centralized servers or peer to peer networksand the like. Different portions of the software programcode/instructions and data can be obtained at different times and indifferent communication sessions or in a same communication session. Thesoftware program code/instructions and data can be obtained in theirentirety prior to the execution of a respective software application bythe computing device. Alternatively, portions of the software programcode/instructions and data can be obtained dynamically, e.g., just intime, when needed for execution. Alternatively, some combination ofthese ways of obtaining the software program code/instructions and datamay occur, e.g., for different applications, components, programs,objects, modules, routines or other sequences of instructions ororganization of sequences of instructions, by way of example. Thus, itis not required that the data and instructions be on a single machinereadable medium in entirety at any particular instance of time.

In general, a tangible machine readable medium includes any tangiblemechanism that provides (i.e., stores) information in a from accessibleby a machine (i.e., a computing device, which may be included, e.g., ina communication device, a network device, a personal digital assistant,a mobile communication device, whether or not able to download and runapplications from the communication network, such as the Internet, e.g.,an I-phone, Blackberry, Droid or the like, a manufacturing tool, or anyother device including a computing device, comprising one or more dataprocessors, etc.

For the purposes of describing and defining the present teachings, it isnoted that the term “substantially” is utilized herein to represent theinherent degree of uncertainty that may be attributed to anyquantitative comparison, value, measurement, or other representation.The term “substantially” is also utilized herein to represent the degreeby which a quantitative representation may vary from a stated referencewithout resulting in a change in the basic function of the subjectmatter at issue.

Although these teachings have been described with respect to variousembodiments, it should be realized these teachings are also capable of awide variety of further and other embodiments within the spirit andscope of the appended claims.

The invention claimed is:
 1. A system for assessing wound, comprising:an image acquisition device configured for capturing an image of a bodypart including a wound area; a segmentation component configured tosegment the image into a number of regions; a region merging componentconfigured to merge regions from the number of regions in order toobtain a smaller number of regions and compensate for over segmentation;a wound boundary determination component configured for extracting aboundary of the wound area; an image color segmentation componentconfigured for performing color segmentation within the boundary of thewound area, wherein the wound area is divided into a plurality ofsegments, each segment being associated with a color indicating ahealing condition of the segment, and a wound evaluation componentconfigured for evaluating the wound area; wherein the wound evaluationcomponent is configured to determine a quantitative healing score as amethod for quantifying a healing status of the wound area; thequantitative healing score being based on the wound area and color inthe wound boundary; wherein the quantitative healing score is a woundevaluation value; the quantitative healing score comprising a firstcomponent determined from a wound area at a predetermined time and aninitial wound area and a second component determined from a ratio ofareas of each color from red, yellow and black; wherein the quantitativehealing score comprises a sum of a product of a first predeterminedweight times the first component and a product of a second predeterminedweight times the second component; wherein the system is used forassessing chronic wounds and ulcers including chronic foot wounds andulcers from type 2 diabetes.
 2. The system of claim 1, wherein thesegmentation component is configured to perform mean shift segmentation.3. The system of claim 1, wherein the wound boundary determinationcomponent is configured to use a trained classifier.
 4. The system ofclaim 1, wherein the image color segmentation component is configured toperform a K-mean color clustering algorithm; and uses a red-yellow-blackevaluation model for evaluation of the color segmentation.
 5. The systemof claim 1, wherein the image acquisition device, the segmentationcomponent, the region merging component, the wound boundarydetermination component, the image color segmentation component and thewound evaluation component are comprised in a handheld portableelectronic device; and wherein the handheld portable electronic devicecomprises: the image acquisition device; the image acquisition devicebeing configured for capturing an image of a body part including a woundarea; one or more processors; and non-transitory computer usable mediahaving computer readable code that, when executed by the one or moreprocessors, causes the one or more processors to: segment the image intoa number of regions; merge regions from the number of regions in orderto obtain a smaller number of regions and compensate for oversegmentation; extract a boundary of the wound area; perform colorsegmentation within the boundary of the wound area, wherein the woundarea is divided into a plurality of segments, each segment beingassociated with a color indicating a healing condition of the segment;and evaluate the wound area; wherein evaluating the wound area comprisesdetermining a quantitative healing score as a method for quantifying ahealing status of the wound area; the quantitative healing score beingdependent on the wound area and on color in the wound boundary.
 6. Thesystem of claim 1, wherein the image acquisition device is comprised ina handheld portable electronic device; and wherein the segmentationcomponent, the region merging component, the wound boundarydetermination component, and the wound evaluation component arecomprised in a computing component; and wherein the handheld portableelectronic device comprises: the image acquisition device; the imageacquisition device being configured for capturing an image of a bodypart including a wound area; one or more first processors; andnon-transitory first computer usable media having computer readable codethat, when executed by the one or more first processors, causes the oneor more processors to: transmit the image to the computing component;and wherein the computing component comprises: one or more secondprocessors; and non-transitory second computer usable media havingcomputer readable code that, when executed by the one or moreprocessors, causes the one or more second processors to: receive theimage from the handheld portable electronic device; segment the imageinto a number of regions; merge regions from the number of regions inorder to obtain a smaller number of regions and compensate for oversegmentation; extract a boundary of the wound area; perform colorsegmentation within the boundary of the wound area, wherein the woundarea is divided into a plurality of segments, each segment beingassociated with a color indicating a healing condition of the segment;and evaluate the wound area; wherein evaluating the wound area comprisesdetermining the quantitative healing score as a method for quantifying ahealing status of the wound area; the quantitative healing score beingbased on the wound area and color in the wound boundary.
 7. The systemof claim 1, wherein the segmentation component, the region mergingcomponent, the wound boundary determination component, and the woundevaluation component are comprised in a computing component; and whereinthe computing component comprises: one or more second processors; andnon-transitory second computer usable media having computer readablecode that, when executed by the one or more processors, causes the oneor more second processors to: receive the image from the imageacquisition device; segment the image into a number of regions; mergeregions from the number of regions in order to obtain a smaller numberof regions and compensate for over segmentation; extract a boundary ofthe wound area; perform color segmentation within the boundary of thewound area, wherein the wound area is divided into a plurality ofsegments, each segment being associated with a color indicating ahealing condition of the segment; and evaluate the wound area; whereinevaluating the wound area comprises determining the quantitative healingscore as a method for quantifying a healing status of the wound area;the quantitative healing score being based on the wound area and colorin the wound boundary.
 8. The system of claim 1 further comprising: animaging component comprising: a first front surface mirror; and a secondfront surface mirror; the second front surface mirror being disposed ata right angle to the first front surface mirror; wherein the imagingcomponent is configured such that an optical path originating at thebody part is reflected at one of the first front surface mirror and thesecond front surface mirror, subsequently reflected at another one ofthe first front surface mirror and the second front surface mirror andterminates at the image acquisition device; wherein the body part ispositioned above the first and second front surface mirrors and awayfrom an axis bisecting the right angle; and wherein the imageacquisition device is positioned above the first and second frontsurface mirrors, away from the axis bisecting the right angle and on anopposite side of the axis bisecting the right angle from the body part.